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ABSTRACT 

We combine strong-lensing masses with SDSS stellar velocity dispersions and HST-ACS effective 
(half-light) radii for 36 lens galaxies from the Sloan Lens ACS (SLACS) Survey to study the mass 
dependence of mass-dynamical structure in early-type galaxies. We find that over a 180-390 km s _1 
range in velocity dispersion, structure is independent of lensing mass to within 5%. This result suggests 
a systematic variation in the total (i.e., luminous plus dark matter) mass-to- light ratio as the origin 
of the tilt of the fundamental plane (FP) scaling relationship between galaxy size, velocity dispersion, 
and surface brightness. We construct the FP of the lens sample, which we find to be consistent with 
the FP of the parent SDSS early-type galaxy population, and present the first observational correlation 
between mass-to- light ratio and residuals about the FP. Finally, we re- formulate the FP in terms of 
surface mass density rather than surface brightness. By removing the complexities of stellar-population 
effects, this mass-plane formulation will facilitate comparison to numerical simulations and possible 
use as a cosmological distance indicator. 
Subject headings: gravitational lensing — galaxies: elliptical 



1. INTRODUCTION 

The fundamental plan e (FP; iDressler et al.l Il987t 
iDiorgovski fc Davis! 1 1987h is a well-known scaling rela- 
tionship between the size R, surface brightness /, and 
velocity dispersion a of elliptical galaxies, expressed in 
the form 

R oc a a I b , (1) 

that indicates an underlying regularity within the pop- 
ulation. Crudely speakin g, the exponents of the FP 
are (a,b) = (1.5,-0.8) fe.g. iBernardi et al.ll2003l ). differ- 
ing significantly from the most nai've constant mass-to- 
light dimensional-analysis expectation of (a, b) = (2, —1). 
This so-called tilt of the FP (relative to the virial plane) 
can be attributed to either a systematic variation in the 
mass-dynamical structure of ellipticals (structural non- 
homology), a systematic variation of the stellar mass-to- 
light ratio, or a systematic variation in the central dark- 
matter fraction with other quantities. This ambiguity 
is difficult to resolve, since neither mass structure (as 
opposed to light structure), nor mass-to-light ratio, nor 
dark-matter frac tion are directly and independen tly ob- 
servable (see e.g. lCiotti. Lanzoni. fc Renzini| [l996L As a 
result, no definite consensus has emerged as to the un- 
derlying explanation of the FP. 
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Strong gravitational lensing can break some of this 
degeneracy observationally, since the total mass within 
the Einstein radius i?Ein is directly observable in grav- 
itational lens galaxies independent of any dynamical 
modeling. The main impediment to this strong lens- 
ing approach has been the lack of a large and ho- 
mogeneous sample of strong lenses for which the tra- 
ditional FP observables R, I, and er can all be re- 
liably measured. This limitation has recently been 
overcome by the Sloan Lens ACS (SLACS) Survey 
(iBolton et all l2006t iTreii et alj 120061 : iKoopmans et all 
120061 : iGavazzi et al.ll2007|: h ereafter Papers I-IV respec- 
tivelv: see also lBolton et ail 2004, 2003), an ongoing sur- 
vey for strong gravitational lens galaxies combining spec- 
troscopic lens-can didate selection fr om the Sloan Digital 
Sky Survey (SDSS: lYorl7et al.l2000D with high-resolution 
follow-up imaging with the Advanced Camera for Sur- 
veys (ACS) aboard the Hubble Space Telescope (HST). 
In this Letter, we investigate the dependence upon mass 
of early-type galaxy mass-dynamical structure by adding 
strong-lensing information to the traditional FP observ- 
ables. We also construct the FP itself for the SLACS 
sample, show that the residuals about the best fitting 
FP are significantly correlated with mass-to-light ratio, 
and present a new formulation of the FP using lensing 
data to replace surface brightness with surface mass den- 
sity. 

2. OBSERVATIONS AND MEASUREMENTS 

The imaging data for all analysis presented in this Let- 
ter were collected with single 420-s exposures through 
the F814W (/-band) filter under HST programs 10174 
(PI: Koopmans) and 10587 (PI: Bolton). As of 2006 
July, 89 systems were successfully observed with ACS- 
WFC by these two Snapshot programs, yielding 44 
strong lenses with redshifts in the range z\ cns ~ 0.1- 
0.4. Full details of the final SLACS Snapshot sam- 
ple will be presented in a forthcoming paper (Bolton 
et al. 2007, in preparation; hereafter B07). To gen- 
erate lens-galaxy-subtracted images suitable for gravi- 
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tational lens modeling and to measure lens-galaxy flux 
within an arbitrary radial aperture, we fit a PSF- 
convolved elliptical b-spline model to the surface bright- 
ness profile of the lensing galaxies (see Paper I). We 
also fit a PSF-convolved elliptical deVaucouleurs sur- 
face brightness model to the image of each lens galaxy, 
to measure intermediate- axis half-light (effective) radii 
(R e ) and brightnesses (I e — 0.5L/irRl) for use in the 
FP analysis. Lastly, we make lensing mass measure- 
ments by fitting sing ular isothermal ellipsoid (SIE; e.g. 
iKormann et al.lll 994f ) mass models to the extended mul- 
tiple images of lensed background galaxies. We cor- 
rect the HST /-b and photometry with the Galacti c 
dust corrections of lSchlegel. Finkbeiner. fe Davis! ff998), 
synthetic fc-corrections to rest-frame l/-band based on 
SDSS spectral templates for each g alaxy, an evolution 
correction of ri l ogTy jdz = —0.4 (jKelson et al.l 120001 : 
iTreu et al.|[200fl : lMoran et al.ll2005f) . and distance moduli 
computed for an (fi M ,ft A ,/i) = (0.3,0.7,0.7) FRW uni- 
verse. Throughout this Letter, the aforementioned FRW 
cosmology is assumed, magnitudes are in the AB system, 
masses and luminosities are in solar units, sizes are ex- 
pressed in kpc, velocities are in krns" 1 , and logarithms 
are base 10. 

In order to achieve a uniform physical aperture with as 
little extrapolation and artificial covariance as possible, 
we correct all lensing mass measurements to an aperture 
of R e /2 using the isothermal mass model, since the me- 
dian ratio of i?Ein to R e within our lens sample is 0.55. 
The isothermal model is preferred by the lens-dynamical 
analysis of Paper III, and is also favored by statisti- 
cal s tudies of lensed quasars (e.g. iRusin fe Kochanekl 
|2005|) and by dynamical analysis of nearby ellipticals 
(e.g. iGerhard et alJl200lh . To assess the sensitivity of 
our results to the isothermal assumption, we derive alter- 
nate aperture-mass corrections using a light-traces-mass 
(LTM) model determined from the b-spline photometry. 
The mean difference between isothermal and LTM aper- 
ture masses within R e /2 is —3%, with an RMS difference 
of ±9%. To compute central mass-to-light ratios (which 
we denote by T e 2), luminosities within R e /2 are mea- 
sured from the b-spline models. We also correct the ob- 
served velocity dispersions from the SDSS fiber radius of 
1"5 to R e j2 (giving g R2 ) using the em pirical relation of 
lJ0rgensen. Franx. fc Kiasrgaardl (|1995l ). These velocity 
corrections are small: 1.4% ± 1.3%. 

We reject 8 lenses based upon (i) prominent spiral 
arms, (ii) an SDSS spectral SNR per pixel of less than 
7, or (iii) significant disagreement between the model R e 
and the radius of an aperture containing one-half of the 
model flux. These cuts leave us with a sample of 36 
early-type strong lens galaxies with absolute evolution- 
corrected magnitudes My in the range —20.8 to —24.0 
and velocity dispersions cr e 2 — 180-390 km s _1 . Quoted 
errors on all fitted parameters are 68% limits from boot- 
strap re-sampling of these 36 lenses. 

3. THE MASS (IN)DEPENDENCE OF STRUCTURE 

Dimensional analysis of the Jeans equation (e.g., 
iBinnev fc Tremaindfl987l ) shows that 

M = ca 2 R/G , (2) 

where R is a scaling radius, a is an average stellar veloc- 
ity dispersion within that radius, M is the mass interior 



to that radius, G is Newton's constant, and c is a dimen- 
sionless structure constant. Note that R and a pertain 
to the luminous (i.e., stellar) component, while M refers 
to the total mass (i.e. stellar and dark): this follows from 
the Jeans equation, and matches well with the quantities 
that we actually measure with photometry, spectroscopy, 
and lens modeling. The value of c will in general de- 
pend upon the anisotropy of the stellar orbits, the over- 
all shape of both the luminous- and dark-matter density 
profiles, and the relative fraction of luminous to dark 
matter in the region of interest. In the SLACS sample, a 
systematic variation of c with galaxy mass will be appar- 
ent in the relationship between the lensing-determined 
total mass -Mi ens within R e /2 on the one hand and the 
dimensional mass variable Md; m = G~ 1 a 2 2 (R e /2) on the 
other, which we parameterize as 

log Mi ens = 8 log M dim + log C . (3) 

A systematic trend of c with mass thus corresponds to 
a value of <5 different than 1. Minimizing the scatter or- 
thogonal to the best-fit relation, we find 6 — 0.986±0.034 
(S = 0.956 ± 0.042) and logc = 0.58 ± 0.16 (logc = 
0.74 ± 0.19) for isothermal (LTM) aperture mass cor- 
rections. This relationship is shown in the left-hand 
panel of Fig. [TJ We see that — whatever the underly- 
ing details — the SLACS lens population is empirically 
consistent with no variation of the structure constant c 
with mass. This result is closely related to the nearly 
one-to-one relation between stellar and isothermal lens- 
model vel ocity disper s ions s hown for the initial SLACS 
sample in ITreu et all ([2006) . By casting this result in 
terms of nearly model-independent aperture masses, the 
current analysis emphasizes the strong evidence for the 
mass independence of early-type galaxy structure in the 
range a e2 =180-390 km s _1 , independent of any assumed 
underlying mass-dynamical model. 

4. PLANES FUNDAMENTAL 

To place the above result in the context of the FP, 
we must first verify that the SLACS sample does indeed 
define the same FP as non-lenses. Paper II showed the 
original SLACS sample of 15 lenses to be consistent with 
the locally determined FP when corrected for luminos- 
ity evolution; the current, larger sample allows us to di- 
rectly constrain the FP of the lenses. We perform this fit 
in the orthogonal sense, defining the best-fitting plane 
as that which minimizes the total squared perpendicu- 
lar distance from all data points in the space spanned 
by logi? e , log/ e , and loga e 2. Defined in this way, the 
coefficients of the FP do not change with the choice of 
dependent variable. Expressing the FP in the form 

log R e = a log cr e2 + b log I e + d , (4) 

we find a = 1.50 ± 0.32, b = -0.78 ± 0.13, and d 

i ■ 1.7. The given errors are the square-root diagonal 
entries of the covariance matrix of a set of bootstrap- 
re-sampled coefficient fits; off-diagonal correlations are 
p ao ~ 0.6, pad — —0.8, and pm — —0.9. The residual 
logarithmic orthogonal scatter about the best-fit plane 
is 0.041 dex. The FP is shown in edge-on projection 
with respect to R e in the center panel of Fig. [T] The 
scaling coefficients (a and b) of the FP defined by the 
SLACS le ns sample are consiste nt with the orthogonal 
FP fits of Bernardi et alj ([2003) for early-type galaxies 
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Fig. 1. — Left: relationship between lensing mass and dimensional mass within R e /2. Dashed line is best fitting linear relation between 
these logarithmic quantities, with slope of 0.986 ± 0.034. Center: effective radius R e as observed and as predicted from the best-fitting FP 
relationship between R e , I e , and cr e 2- Right: mass-to-light ratio T e 2 in the V-band versus residuals about the FP relationship. 



from the SDSS (with somewhat different selection and 
conventions) . 

As an aside, the addition of lensing data allows us to 
make the first direct observation of a correlation between 
the residual scatter about the FP and the total mass-to- 
light ratio T e2 as determined from lensing, shown in the 
right-hand panel of Fig. [1] The correlation has a linear 
(Pearson) coefficient of —0.58, giving a formal signifi- 
cance of 99.98% for a sample size of 36. The significance 
of the correlation is similar under the LTM assumption 
(i.e. no correction of T from i?Ein to R e /2). Furthermore, 
the FP residuals do not have any significant correlation 
with mass or with luminosity separately, nor any cor- 
relation with redshift or with the ratio R^i n /R e . We 
are thus drawn to the conclusion that this correlation 
is intrinsic, and not introduced through a correlation of 
measurement errors. Though they did not observe it 
directly, such a correlation o f FP residuals wit h mass-to- 
light ratio was deduced by Fa ber et ai1ll987l to be the 
most likely ex planation for the intrinsic thickness of the 
FP. (Also see lJ0rgensen. Franx. fc Kiaergaardlll996l for 
an extensive observational analysis of FP residual corre- 
lations and intrinsic thickness.) 

We now use the strong-lensing data to formulate an 
analogous plane in mass space, which we will refer to as 
the mass plane (MP). Following the traditional form of 
the FP, we define the MP by replacing surface brightness 
I e with surface mass density within i? e /2, denoted by E e2 
(inM©kpc~ 2 ): 

log R e = a m log f7 e2 + b m log £ e2 + d m . (5) 
Note that in going from I e to E e2 , the sensitivity to 
luminosity evolution with redshift is removed. Fitting 
as for the FP above, the MP coefficients are a m = 
1.77 ± 0.14, b m = -1.16 ± 0.09, and d m = 7.8 ± 1.0. 
Off-diagonal correlations are p a b,m — 0.1, p a d,m — —0.4, 
and Pbd,m — —0.9. The RMS orthogonal logarithmic 
scatter about the MP is 0.026 dex. Fitting for the MP 
with mass densities computed using the LTM assumption 
rather than the isothermal model gives a m = 1.86 ±0.17, 
b m = —0.93 ± 0.09, and d m = 5.4 ± 0.9, with compara- 
ble parameter correlations and an orthogonal logarithmic 
scatter of 0.030 dex. The MP thus appears to be tighter 
than the FP in the sense of having smaller residual scat- 
ter. This scatter is consistent with the ~7% fractional 
velocity-dispersion errors in the current data. 

The significant difference between the FP and the MP 
is most easily seen by considering the parameterization 



of the FP tilt i n terms of a relationship of the form 
L oc MZ m (e.g. iDrgssler et al l [Ml; ICiotti et"aLl ll996L 
ITruiillo. Burkert. fc Belll200l iTreu et al.ll2005bf ) . While 
this form does not exactly reproduce the best-fitting FP 
exponents of the current study, it is a convenient phys- 
ically motivated approximation. Fitting in log-space by 
minimizing the orthogonal scatter, we find r\ = 0.82±0.05 
for the SLACS sample. This translates into FP coeffi- 
cients of a — 1.39 and b = —0.85 in the form of Eq. |U 
Comparing this with the result that 6 = 0.986 ± 0.034 
from § [3] — which amounts to an analogous parameteriza- 
tion of the MP — we see that in going from the FP to the 
MP, the tilt relative to the virial relation is essentially 
eliminated. 

5. DISCUSSION 

Previous studies have co nstructed the FP of 

gravitational lens ga l axies (fKochanek et alj 120001 : 
Ivan de Yen et all l2003t iRusin fc Kochanekl l2005h . but 
have taken their velocity dispersions from the strong 
lensing data rather than from stellar dynamics. By 
combining independent measurements of mass and 
velocity, we have found direct evidence that the total 
mass-dynamical structure of early-type galaxies is 
independent of mass. What, then, does this result 
tell us about the underlying explanation for the tilt 
of the FP? The straightforward binary answer is that 
the FP must be due to a systematic variation in the 
total central mass-to-light ratio — either through the 
stellar population or through the central dark-matter 
fraction — and not due to any systematic variation of 
mass-dynamical structure with mass. If mass-dynamical 
structure effects were responsible for the FP, they should 
also manifest as a systematic variation with mass of the 
structure constant c of § [3] and a similar tilt in the MP 
as in the FP, whereas we see neither of these two effects. 
This conclusion is most no tably consistent with those o f 
iPadmanabhan et all (|2004h and lCappellari et all (|2006h : 
the agreement is reassuring given the diversity of the 
methods employed between these studies and ours. 

Our conclusions are at variance with those of 
iTrujillo et all (|2004h . who point to non-homology in the 
luminous component of early-type galaxies (specifically, 
the systematic variation in Sersic index with galaxy lu- 
minosity) as being largely responsible for the FP. A 
pl ausible explanation i s that , motivated by t he res ults 
of iRomanowskv et al l (|2003D . ITruiillo et ail (|2004D do 



4 



Bolton et al. 



not allow for the presence of dark matter. Other lines 
of evidence suggest a "bulge-halo conspiracy" in early- 
type galaxies, whereby the combined luminous and dark- 
matter profiles generate a more uniforml y flat rotation 
curve than either component alone (e.g. I Gerhard et al.l 
|2001[ ) . Such an effect could offset non- homology in the lu- 
minous component, leading us back to systematic mass- 
to-light variation as the explanation for the FP. 

As mentioned above, our observations do not distin- 
guish between stellar mass and dark mass: the two physi- 
cally very distinct possibilities of a variation in the stellar 
mass-to-light ratio and a variation in the central dark- 
matter fraction are (in principle) equally allowed. Pre- 
vious studies have shown that variations in stellar mass- 
to-light must be responsible for some (but not all) of the 
tilt of the FP, as seen through the systema t ic dec rease 
in tilt at redder wave-bands dPahre et alj 119981). de - 
tailed stellar population modeling ( Thomas et al.l l2005). 
and the galaxy-mass-depende nt average rate o f lumi- 
nosity evolution with redshift (|Treu et al.ll2005aH bT). At 
this stage, the main diagnostic available to us to quan- 
tify stellar populations effects in the SLACS sample is 
the strength of the 4000A break in lens-galaxy contin- 
uum flux £)„(4000), w hich should co r relate with stel- 
lar mass-to-l i ght ratio (jBruzuall 119831: [Hamilton! 119851 : 
iBalogh et all 119991 : IKauffmann et all 120031) . No signif- 
icant correlation with T e2 is found, although the sen- 
sitivity of Z?„(4000) may not be sufficient to detect 
the small age differences (Ai/f ~ 2/10) expected for 
early-type galaxies i n this range of velocity dispersion 
(|Thomas e t al. 2005). Taken at face value, the result 
suggests that most of the variation in T e 2 is due to dark 
matter content. We will quantify this statement in future 
work with the aid of more sophisticated stellar popula- 
tion diagnostics. 

In the process of investigating the mass dependence 
of early-type galaxy structure and the genesis of the 
FP, we have made the first observational correlation of 
the thickness of the FP with mass-to-light ratio, and 
we have introduced an analogous plane in mass space 
that is both tighter and less tilted than the FP relative 
to the virial scaling relation. This observational mass 
plane is a fundamental dynamical scaling relation for 
early- type galaxies; though it requires a measurement 
of galaxy size, it makes no reference to galaxy lumi- 
nosity, and therefore it can be tracked across cosmic 
time without regard to stellar evolution. The MP thus 
offers great promise as a cosmological standard ruler, 
particularly if the relation can be calibrated locally by 



masses from detailed dynamical models of nearby ellip- 
tical galax ies. (Though structural evolution m ay still 
occur: e.g. Ivan der Marel &: van Dokkumll2^0falrbl .) By 
removing the need for stellar-population modeling, the 
MP is also more suitable for comparison with the results 
of numerical simulations. 

This Letter has not addressed possible biases intro- 
duced through our gravitational-lens selection procedure. 
We defer a full discussion of selection effects to B07, but 
note that Papers I and II showed the original SLACS 
sample to be statistically consistent with the parent sam- 
ple of non-lens SDSS galaxies from which they were se- 
lected. 

The analysis presented here has been made possible 
only recently by the SLACS Survey gravitational lens 
sample. Existing and forthcoming deep HST imaging 
will permit more accurate mass, luminosity, and color 
measurements, while ongoing spectroscopy with the VLT 
and Keck telescopes will provide more robust velocity 
dispersion measurements and stellar-population diagnos- 
tics. These data will soon enable a greater precision in 
the type of analysis described in this Letter. 
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